Local Properties of the Potential Energy Landscape of a Model Glass: Understanding 

the Low Temperature Anomalies 
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Though the existence of two-level systems (TLS) is widely accepted to explain low temperature 
anomalies in the sound absorption, heat capacity, thermal conductivity and other quantities, an 
exact description of their microscopic nature is still lacking. We performed computer simulations for 
a binary Lennard- Jones system, using a newly developed algorithm to locate double-well potentials 
(DWP) and thus two-level systems on a systematic basis. We show that the intrinsic limitations 
of computer simulations like finite time and finite size problems do not hamper this analysis. We 
discuss how the DWP are embedded in the total potential energy landscape. It turns out that 
most DWP are connected to the dynamics of the smaller particles and that these DWP are rather 
localized. However, DWP related to the larger particles are more collective. 



INTRODUCTION 



It is well known that most kinds of disordered solids 
show anomalous behavior at very low temperatures as 
compared to their crystalline counterparts. Many of the 
observed features can be explained by the Standard Tun- 
nehng Model (STM) 0, and its generalization which 
is the Soft-Potential Model The basic idea 

of the STM is the existence of a broad distribution of 
Two Level Systems (TLS). The TLS can couple to strain 
and electric fields and therefore influence quantities like 
the heat capacity, thermal conductivity, sound absorp- 
tion and dielectric response, see for a review. The 
STM predicts a linear dependence of the heat capac- 
ity on temperature and a quadratic dependence of the 
thermal conductivity on temperature. The STM gives a 
good general agreement with experimental results down 
to temperatures around 100 mK. 

So far it has not been possible to derive a theory of the 
glass transition or of the low-temperature anomalies from 
first principles, i.e. from the Hamiltonian of the glassy 
system. This means, that, except for the few cases where 
possible TLS have been successfully identified by means 
of computer simulations |^ 0| , the STM and the recent 
developments, mentioned above, are almost purely phe- 
nomenological. In particular researchers wanted to ex- 
plain why the nearly constant ratio of the density of TLS 
and their coupling to phonons, if compared for very differ- 
ent glasses, is so similar [13. Already one decade ago this 
has been interpreted by Yu and Leggett 0,^3 ^'^^ Cop- 
persmith Jjj as an indication that the observed TLS are 
highly collective excitations of many underlying micro- 
scopic TLS, resulting from the interaction among TLS. 
In the meantime, however, it has been shown that the 
interaction is only relevant in the mK regime [13. ITsIITgII . 
An alternative scenario has been proposed by Lubchenko 
and Wolynes JTJ. They consider the glass as a mosaic 
of frustrated domain walls, separating individual cells. 
In their model the collective tunneling process finally in- 



volves O(IO^) molecules which only move a fraction d^/a 
of a nearest-neighbor distance {{dL/aY k, 0.01). Unfor- 
tunately, the object of these theories, namely the TLS, 
are still somewhat obscure because experimentally it is 
very difhcult to characterize their microscopic nature. 

Formally, a TLS corresponds to a pair of local min- 
ima, or a double- well potential (DWP), on the potential 
energy landscape (PEL). The minima need to have an en- 
ergy difference less than fc^T and a small distance in con- 
figuration space because otherwise no tunneling would 
occur. One may wonder whether computer simulations 
might help to elucidate the relevant properties of TLS 
and thus to prove their existence. Computer simulations 
are strongly limited in several directions: (a) Typically 
rather small systems have to be used to analyze the PEL. 
This may give rise to significant finite size effects; (b) 
Due to finite simulation times the cooling process is ex- 
tremely fast so that the resulting glassy structure may 
be vastly different as compared to the experimental situ- 
ation. This may strongly influence the properties of the 
TLS obtained by computer simulations; (c) Due to possi- 
ble imperfections of the search algorithm to identify TLS 
one may possibly miss a significant fraction of TLS. 

In this contribution we will show that these problems 
are not relevant for the problem of locating TLS and 
that it is indeed possible to obtain detailed and unbiased 
information about the nature of TLS. A particular chal- 
lenge is the systematic search of TLS. Most of the early 
computer simulations in this field 0, 0, [23, [^ [IJ 
did not attempt to systematically find DWP. In this 
work we present computer simulations on a model glass 
former (binary Lennard-Jones) to systematically iden- 
tify the TLS. A first step in this direction has been 
already published a decade ago [2^ [3j [IS- At that 
time, however, it was not possible to exclude that any of 
the above-mentioned problems might hamper the anal- 
ysis. With improved algorithms and faster computers 
this has become possible nowadays. Furthermore the 
TLS are related to the properties of the glass transi- 
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tion. Qualitatively, one may say that TLS probe the 
PEL on a very local scale whereas for the understanding 
of the glass transition much larger regions of the PEL 
are relevant. Here we would like to mention that in re- 
cent years computer simulations succeeded in extracting 
many important features of the PEL of supercooled liq- 
uids m El El El US El HI mm. 

II. TECHNICAL 

A. Computational details 

As a model glass former we chose a binary mixture 
Lennard-Jones system with 80% A-particles and 20% B- 
particles (BMLJ) [H El IP, "stl . It is supposed to rep- 
resent Nickel-Phosphorous (80% 62Ni; 20% ^ip) jgj but 
with a 20% higher particle density, this system was first 
used by Kob and Anderson. The used potential is of the 
type 

V^a/s = 4 • ea/3[K;3/r-)i' - {^c^p/rf] + (a + 6 • r), (1) 

with aAB = O-SaAA, ctbb = O.SSaAA, ^ab = l-Se^iA, 
esB = O.SeAA, TTT-B — O.Stoa- Periodic boundary con- 
ditions were used and the linear function a + b ■ r was 
added to ensure continuous energies and forces at the 
cutoff Tc = 1.8. The units of length, mass and en- 
ergy are aAA, "m-A, ^aa, the time step within these units 
was set to 0.01. The simulation cell was a cube with 
a fixed edge length according to the number of parti- 
cles and an exact particle density of D = 1.2. Molecu- 
lar dynamics (MD) simulations, using the velocity Ver- 
let algorithm, have been used to generate independent 
configurations and as part of the DWP location algo- 
rithm. For the case of Nickel-Phosphorous the energy 
unit corresponds to 933.9 K and ctaa is 2.2 A. We an- 
alyzed = 65, 2 • 65, 130, 2 • 130, 195, 260 particle 
systems. The 2 • 65 and the 2 • 130 systems denote two 
non interacting systems sharing one simulation box. 

It is known for the same system that the dynamics 
above the mode-coupling temperature Tc as well as struc- 
tural properties like the pair correlation function arc ba- 
sically independent of system size for A^ > 65 _29, 38j . 
Therefore one may hope that even for the very small sys- 
tems finite size effects for the TLS can be neglected. It 
will turn out that this is indeed the case. 



B. Systematic Location of TLS 

We have developed a new algorithm for a systematic 
search of TLS. Formally, the problem is to identify two 
nearby local minima on the potential energy landscape 
(PEL) of the system. In the first step a set of equi- 
librium configurations is generated via MD simulations 
at constant temperature T^quU- Configurations at differ- 
ent times are taken and used as starting configurations 



for subsequent minimization via the Polak-Ribiere conju- 
gate gradient algorithm, yielding a corresponding set of 
local minimum energy structures, denoted inherent struc- 
tures . In the second step the goal is to locate nearby 
minima for these inherent structures. To search for these 
nearby minima we proceed as follows: (1) Reset all parti- 
cle velocities with random numbers according to a normal 
distribution at a fixed temperature Tsearch- (2) Perform 
a fixed number Ngearch of MD steps in the NVE ensem- 
ble. (3) Minimize the resulting structure. (4) Accept the 
new minimum if the distance d and the asymmetry A 
between both minima in configuration space (see below 
for an exact definition) fulfill ~ dmin < d < dmax and 
« Amin < c^A < A„jaa;, respectively. We introduced 
very small minimum cutoffs, because minima cannot be 
distinguished below numerical precision. For the choice 
of the cutoff- values, see below. (5) Repeat this procedure 
M times for different initial random numbers. The val- 
ues of Tsearch and N search are choscu to give the largest 
number of nearby minima per starting inherent structure. 

In the final step the transition states for the generated 
pairs of minima are located by a modified version of the 
non-local ridge method '38]. This algorithm is very ro- 
bust to identify first-order saddles between pairs of min- 
ima. In case that no such saddle exists the new minimum 
is dismissed. In particular this may occur if the path to 
the second minimum involves transitions of two (or more) 
independent TLS during the Nsearch MD-steps. This is 
the standard situation for very large systems. As a result 
we obtain for every starting minimum a set of minima 
characterized by asymmetry A, distance d and barrier 
height V . We explicitly checked that all minima have 
and all transitions states 1 negative eigenfrequency. 

This algorithm is superior to the TLS-search algo- 
rithm, used in previous works [HE^ because no indirect 
assumptions are made with respect to the number of par- 
ticles participating in a TLS-transition or with respect to 
the distance d. 

Beyond the direct euclidian distance d between config- 
urations (or correspondingly , the distance di moved by 
particle i), we use mass weighted distances dmw between 
two configurations 

N 

Cjn,ri) = ^(dL + + dl) ■ ^. (2) 

i 

Furthermore one can define the mass weighted reaction 
path approximation between two minima via 

dmwrp — dyyiyj{T\, Tiri^dyissidie^ ~\- drnwiXtrans. states ^2) • (3) 

Where fi,f2T'r'tra7is. state are the particle positions of the 
two minima and the transition state, drp is defined sim- 
ilarly as dmwrp 1 but without mass- weighting. The mass- 
weighted reaction path is introduced because it enters the 
WKB-term to calculate the tunneling matrix element. 
For comparison with literature it is necessary to have also 
direct euclidian distances. In the subsequent analysis we 
have attempted to employ for every type of analysis the 
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appropriate distance. In any event, due to the similarity 
of all definitions, none of our results would change on a 
qualitative level if different distances had been used. 



C. Double-well potentials vs. TLS 

As the ultimate goal one is looking for pairs of min- 
ima, i.e. double- well potentials (DWP), in the high- 
dimensional potential energy landscape with A < IK, 
corresponding to 0.001 in LJ-units. They are expected 
to be relevant for the low-temperature anomalies around 
IK and thus act as TLS. Qualitatively, one would ex- 
pect that many more DWP exist with asymmetries of 
the order of the glass transition temperature rather than 
IK. Restricting the search to DWP with such a small 
asymmetry is somewhat problematic because TLS are a 
very rare species and thus it is very difficult to find them 
numerically. Therefore we use as a standard choice a rel- 
atively large asymmetry of Amax — 0.5. From the result- 
ing data set (6522 DWP for N = 65, 2911 for TV = 130 
and 428 for N = 260) we may easily analyze subsets 
of almost symmetric DWP which are relevant for the 
low-temperature properties. As shown in previous work 
also the range of very asymmetric DWP contains impor- 
tant information about the properties of TLS. In partic- 
ular the choice A^ax = 20 K (0.02 in LJ-units) is small 
enough as compared to the glass transition temperature 
but large enough so that we get sufficient statistics. We 

0.8 for reasons mentioned be- 



mwrp^max 



always choose d 
low. For the selection of the starting minima we choose 
Tequii = 0.5 which is slightly above the critical mode- 
coupHng temperature of Tc= 0.45 [sslls^. 



III. RESULTS 

A. Finite size effects and completeness of the 
search 

From previous work it is known that the different pa- 
rameters A, d, and V are strongly correlated. For ex- 
ample one finds in agreement with intuition that DWP 
with small d typically have a rather small asymmetry. 
Thus one would expect that most DWP which may act 
as TLS at low temperature are restricted to some re- 
gion of d- values. Therefore we have checked whether the 
DWP with A < 0.02 indeed display this restriction in 
d. For = 65 we have found 301 DWP in this range. 
The result is shown in Fig. ^ It turns out that the d- 
distribution shows a peak around c? = 0.3 and strongly 
decreases for larger d. For A^ — 130 we get very similar 
results. Therefore our choice dmax — 0.8 can be justified. 
In a next step we elucidate the quality of our search algo- 
rithm and check whether our search of TLS is complete 
within our specified parameter range. This property is 
essential to estimate the absolute number of TLS from 
our simulations. For this purpose we have analyzed how 
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Figure 1: The histogram of the distribution of distances for 
all DWP with A < 0.02 (20K for NiP). 
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Figure 2: The histogram of how often a nearby minimum was 
found for the 65 particle system after M = 400 attempts per 
starting minimum. The triangles correspond to data limited 
to participation ratios above average (1/3 of all found). 



often the different minima around a starting minimum, 
i.e. how often DWP, are found during the M = 400 at- 
tempts. In case that DWP are only found once or twice it 
is very likely that many minima are overseen. In contrast, 
if all DWP are found quite frequently it is unlikely that 
many DWP are overseen. This conclusion is based on 
the assumption that properties of DWP are continuous, 
i.e. there is not a second class of invisible DWP, strictly 
separated from the DWP found in our simulations. For 
A^ = 65 particles the distribution of TLS counts is shown 
in Fig. 121 for both A^ax = 0.5 and A^ax = 0.02. Al- 
ready for the first choice it turns out that the maximum 
of the distribution has its maximum at a count above 
15 and the likeliness for lower counts decreases rapidly. 
Thus most minima are indeed found. For the DWP with 
A < 0.02 this effect is even more pronounced. This com- 
parison also shows that DWP with smaller asymmetry 
are found more easily than DWP with larger asymmetry. 
Furthermore we analyzed whether for the system size 
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of iV = 65 finite size effects are present. As already 
mentioned above, finite size effects are absent for the 
dynamics above Tc- It turns out tliat all properties of 
DWP (partly presented in this work, partly in a sub- 
sequent publication) ware identical in the observed N- 
range between 65 and 260 within statistical noise. There 
is, however, one exception. Within the parameter range, 
specified above, we found on average 0.65 DWP per start- 
ing minimum for TV = 65 and 0.94 DWP for N = 130. 
We used M = 400 and M = 800, respectively. Naively 
one would have expected twice the number of DWP for 
N — 130. This means that on the per-particle basis the 
number of TLS is roughly 30% too smaU for N = 130 
as compared to = 65. This effect is even more pro- 
nounced for N = 260. Two reasons are possible: either 
we have found a significant finite size effect (which may 
appear surprising because all properties of the TLS are 
identical) or this effect may be caused by the algorithm. 

To test the second possibility we performed simula- 
tions with two non interacting systems A and B sharing 
a box. For determining minima the relative positions 
of the two independent systems were preserved and the 
energies added. In the combined system the absolute 
number of TLS is by construction twice as high as in the 
elementary system. Interestingly, using the same simu- 
lation parameters as above we also observed a relative 
decrease of circa 30% when analyzing the combined sys- 
tem 2-65. Thus it is very likely that the apparent finite 
size effect, described above, is due to the algorithm. This 
effect can be rationalized. If a minimum, corresponding 
to a TLS, is found in configuration A with a probabil- 
ity of 90% per attempt and in configuration B another 
TLS is found with a probability of 10% per attempt, they 
change to 89% and 1%, respectively, when both config- 
urations are considered together and if combined transi- 
tions are ignored (see above). This effect gives rise to a 
sharp decrease in probability to detect elementary tran- 
sitions for larger systems. Therefore the best procedure 
is to analyze rather small systems which, however, are 
large enough to be void of relevant finite size effects. In 
analogy to our results in [s^ the system size of A^ = 65 
seems to be a very good compromise and, according to 
Fig. 121 allows one to obtain the absolute number of TLS. 

To proceed we counted the number of DWP with 
A < 0.02. It turns out that one DWP exists per 1000 
particles. For even smaller asymmetry A < IK (« 0.001) 
we observe one DWP, i.e. TLS, per 12000 particles. Us- 
ing the density of NiP one ends up with 6 • 10''^ J~^to~^ 
TLS. Actually, the parameter P, which can be deter- 
mined experimentally and gives an experimental num- 
ber of TLS, is roughly smaller by one decade because 
the contribution of the individual TLS are weighted by 
a factor Ao/(Ao -I- A^) < 1 where Aq is the tunnel- 
ing matrix element. More specifically we obtain P k, 
1.3 • lO'^V-^m-^liJ which is close to the value of P « 
1.5 • lO'^'^J^^m"^, obtained in [IJ. Most of the remain- 
ing difference with, e.g. P = 8 • 10'*'* J~^m~'^ for silicate 
can be directly explained by the fact that most glasses 



have larger elementary units like the Si04 tetrahedra for 
silicate. This means that the number of DWP per vol- 
ume is even smaller, yielding an estimate very close to 
the experimental value of silicate. 



B. Embedding of TLS in the potential energy 
landscape 

Now we discuss how the TLS are embedded into the 
overall PEL. As already discussed above, computer sim- 
ulations suffer from finite simulation times. As a conse- 
quence the starting inherent structures will have a rel- 
atively high energy as compared to inherent structures 
the system would have reached if equilibration at lower 
temperature TequU close to the calorimetric glass tran- 
sition were possible. In order to elucidate this aspect 
closer we would briefly like to summarize the previous 
results for this LJ-system [s^- (1) Inherent structure 
energies range between -304 and -287. The low-energy 
cutoff, however, is not due to the bottom of the PEL. 
Rather it indicates that the number of minima with even 
lower energy is so small that they were not detected dur- 
ing the simulations. (2) More strictly, the distribution of 
energies follows the left part of a Gaussian distribution. 
This implies that the number of states with higher en- 
ergy is always exponentially larger than the number of 
states with somewhat lower energy. (3) It can be esti- 
mated that the bottom of the potential energy landscape 
is around -306. This number can be estimated from de- 
termination of the total configurational entropy [i^ l43l| . 
Thus even for very small cooling rates energies would be 
larger than -306. (4) The dynamics around and below Tc 
can be interpreted as jumps between different traps, de- 
noted metabasins. They contain a number of individual 
inherent structures. 

As our simulations cover a broad range of minimum 
energies we checked to which degree the properties of 
the DWP depend on the energy of the starting inherent 
structure. Fig.Oshows the dependence of the three DWP 
parameters on this energy. Interestingly, the dependen- 
cies are rather small. In particular there is no indication 
to believe that DWP in the relevant low-temperature en- 
ergy range between -306 and -302 are very different. This 
is in agreement with the observed independence of the in- 
herent structure dynamics observed by Vogel et al. |33j |. 
Of course, on a strict quantitative level minor variations 
should be taken into account. 

Furthermore we analyzed whether the number of DWP 
per starting configuration depends on its potential en- 
ergy, see Fig. 0J In order to work with a larger data 
set we have included asymmetries up to 0.8. Further- 
more we have analyzed the DWP for different subsets 
with respect to the distance dmwrp- It turns out that the 
dependence on energy is very weak. Even in the most 
extreme case {dmwpr < 0.8) the number of DWP changes 
only by a factor of 2 when comparing different energies. 
This small variation nearly vanishes when restricting our- 
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Figure 3: The dependence of the DWP parameters on the 
energy of the starting minimum. The lines are a guide to the 
eye. 
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Figure 4: The dependence of the number of DWP per starting 
minimum on the energy of the starting minimum. 



selves to DWP with small distances. In summary, we may 
conclude that on a local scale the properties of the po- 
tential energy landscape do not depend on the height in 
the landscape. Evidently this property must break down 
when taking into account pairs of minima with larger 
spatial separation. Because the number of minima for, 
e.g., -295 is by a factor of 10^ larger than the number 




Figure 5: One dimensional sketch of the potential energy land- 
scape. The asymmetry Adwp and the barrier height for the 
DWP vary between and 2 Tc- The barrier height for the 
traps Vrrap on the other hand varies between 1 and 20 Tc- 



of minima for —300 it is evident that for larger d many 
more pairs of minima with similar energy can be found 
for -295 as compared to -300. 

In Fig. El we attempt to sketch the PEL of the LJ- 
system. The most prominent features are the individual 
traps. In contrast, the DWP analyzed in this work cor- 
respond to the little wiggles within the traps. Typical 
energy scales have been given in the figure caption show- 
ing the separation of transitions between adjacent inher- 
ent structures, already relevant at low temperatures, and 
between traps, relevant above the glass transition. 



C. Microscopic nature of DWP 

Furthermore we have studied the microscopic nature of 
the DWP in great detail. As one of the most elementary 
questions one may ask how many particle are involved 
in the translation between the two minima of a DWP. 
To establish an accurate picture of the microscopic na- 
ture of the TLS we analyzed the number and types of 
particles participating in the transition from one mini- 
mum to the other. As a matter of fact all particles move 
when going from one minimum to the other. In Fig. we 
show the distribution of distances moved by the different 
particles during the transition between both minima of a 
DWP. The distances are sorted according to their actual 
size. The curvature for the last particles is caused by 
the minimization of the distance of the two minima in a 
finite system. The curves for different systems show only 
small deviations for the particles which move a larger dis- 
tance. It can be seen that the first particle, which moves 
the biggest distance, moves much further than the second 
particle if it is a B-particle (small), this is less pronounced 
if it is an A-particle. As the particle which moves the 
biggest distance can usually clearly distinguished from 
the other particles we denote it as central. With this def- 
inition we find that in over 90% of all DWP the central 
particle is a B-particle, although the concentration of B- 
particles is only 20%. DWP with central B-particles thus 
dominate DWP with central A-particle. 

From these data one may define a participation ratio 
to estimate the number of particles involved in the tran- 
sition between two minima. In case that n particles are 
moving some fixed distance and the other particles do 
not participate at all one would obtain a participation 
ratio of n. In the present case, of course, one has a broad 
distribution and thus different definitions of the partic- 
ipation ratio yield different values, as shown in Tab. 
The value for the total system varies between 2.1 and 
9.5. Therefore it is more informative to consider the to- 
tal distribution of translations, as shown in Fig. 

Finally, one may ask the question whether the partici- 
pation ratio depends on the properties of the DWP. Most 
naturally, one might think of a relation to the distance be- 
tween both minima because the distance directly appears 
in all definitions of the participation ratio. In Fig. [3 the 
dependence of the participation ratio (second definition) 
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particle 
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(^d 1 dmax ) 


^ ^ di j dmax ) 




average 


2.14 


3.29 


9.50 


7.65 


B 


2.04 


3.01 


8.97 


6.95 


A 


4.95 


6.72 


16.11 


16.43 



Table I: Partition ratios for different definitions and for A and 
B-particles as central particles. No mass-weighting has been 
used for the distances. 
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Figure 6: The contribution of each particle to the total trans- 
lational motion in a DWP. The DWP with an A-particle as 
central particle (bottom) and a B-particle as central particle 
(top) show different behavior. 



on the distance is shown. As expected for DWP with 
larger distances more particles contribute to the trans- 
lation. Actually, in previous work ("Z?] P. 483) we have 
shown for a very similar LJ-system that the slope might 
decrease for larger distance, still the participation ratio 
increases with distance. 
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Figure 7: The participation ratio does significantly dependent 
on the distance of the configurations. 



D. Spatial distribution of DWP 

In the STM it is assumed that the TLS are randomly 
distributed in space. To check this hypothesis we first 
calculated the probability to find two DWP which have 
the same central particle. We find that, for both N — 65 
and N = 130 an enhanced probability that a single par- 
ticle is central for two different DWP. To elucidate this 
aspect further we made a simplified analysis for the 65 
particle system. Instead of considering all pairs of central 
particles we considered only B-B-pairs and only those 
pairs, which originate from configurations with exactly 
two observed DWP. For independent DWP one would 
expect that in 1/13 — 7.7% of all cases the central parti- 
cle is identical. Rather we found (17 ± 2)%. Thus there 
is indeed a significantly increased probability that two 
DWP are spatially correlated, which can also be inter- 
preted as an increased probability for triple well poten- 
tials. We find that when omitting DWP with the same 
central particle the other DWP are randomly distributed 
in the system. This was checked by computing the av- 
erage distance of the remaining central particles in cases 
where more than 1 DWP is present in the system. It 
turns out that this distance is within statistical uncer- 
tainty identical to the distance of randomly chosen par- 
ticles. The results for the 130 particles system show the 
same behavior. 



IV. DISCUSSION AND SUMMARY 

We have presented a new reliable algorithm to system- 
atically locate DWP in a model glass former. It turned 
out that the intrinsic limitations of computer simulations 
do not hamper the quality of our results. Thus we can 
indeed get information about TLS, relevant for under- 
standing the low-temperature anomalies. 

It turns out that the number of TLS, directly obtained 
from our data, is compatible with the number of TLS ob- 
served experimentally. This conclusion had been already 
drawn from our previous work. This time, however, we 
can exclude that possible systematic artifacts hamper our 
analysis. Thus it is likely that TLS, responsible for the 
low-temperature anomalies, result from elementary two- 
state systems rather than from collective excitations. 

As a by-product of our simulation we realized that the 
occurrence of DWP is not randomly distributed in space. 
Rather there is a strong tendency that two DWP are 
located at the same central particle, i.e. form a triple 
well potential. This suggests that some structural fea- 
tures of the glass favor the formation of DWP. A possible 
candidate is the coordination sphere around the central 
particle. Whether or not these spatial correlations are 
still significant for DWP with asymmetries in the Kelvin 
range would require a much larger set of DWP than can 
be found by current simulations. 

The properties of the TLS reflect local properties of 
the PEL. In connection with our previous work on global 
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features of the PEL it is possible to have a view on the 
PEL, encompassing the transport dynamics above and 
the local dynamics below the glass transition. It becomes 
evident, that the TLS are part of the individual traps in 
which the PEL can be decomposed. 

We would like to mention two discrepancies with other 
work. One is the magnitude of the absolute movement of 
a central B or A-particle, which is found to be quite large: 
r-p — ^-^"^ ^"a rp — 0-06 in LJ units. These values 
are about an order of magnitude larger than the esti- 
mation given by Lubchenko and Wolyness [T^ in their 
frustrated domain wall model. The other discrepancy 
concerns the participation ratios. The observed partici- 
pation ratios, obtained from averaging over all DWP, are 
much lower than those observed by other groups inves- 
tigating the PEL, namely by Oligschleger and Schober 
in soft s phe re glasses and LJ-systems 0, [s^ or Vogel 
et. al. |33l| in the same system as used in this work. 
In [sil l32l 33] the authors have analyzed the transitions 
between adjacent minima as resulting from a molecular 
dynamics trajectory at a given temperature and found of 
the order of 20 particles. Comparing this with our aver- 
age value of around 3 this seems to be a major difference, 
as already stated by Oligschleger and Schober [sj] in re- 
lation to our earlier work where similar small values 
have been reported. As shown in Fig. |21 this difference 
is not due to the fact that our algorithm is not able to 
identify DWP with large participation ratios. For a closer 
discussion of this discrepancy one has to take care of the 
actual definition, used to characterize the participation 
ratio. In the reported work the two latter definitions in 
Tab. m have been used. Thus the reported values have 
to be compared to our participation ratios of about 7-10. 
There is, however, still a remaining difference of a factor 
of 2-3. 

This difference can be rationalized by the different 
methods of locating inherent structures. In the present 
work we have attempted to localize all inherent structures 
within a certain distance to the original minimum. This 
approach was motivated by the observation that DWP 
with small asymmetries typically correspond to nearby 
minima (see Fig.l) so that DWP with large distances be- 
tween the minima are typically irrelevant for the under- 
standing of the low-temperature anomalies. As a direct 
consequence the typical distances between the minima 
are significantly smaller found by the present method as 
compared to the results reported in [s^ which 
more allude to the properties of the glass transition. 
This can be quantified via the configurational distance 
= J2iidi,x + di^y + ^)°'^- Comparing our value with 
the value reported in [sj we find a difference of a factor 
of three. 

Following the results, reported in Fig.Q this difference 



directly translates into differences of the participation ra- 
tio. Extrapolating the data shown in Fig. [7|to distances 
around unity one ends up with an increase of the partic- 
ipation ratio by a factor of 2-3 which is exactly the fac- 
tor, which was missing above. Thus our present results 
are fully consistent with the results discussed in previous 
work. We just mention in passing that the situation may 
be even more complicated since the DWP, obtained by 
Schober et al. and Vogel et al., have been obtained from 
MD trajectories. Thus there exists an implicit weight- 
ing by the probability to find these DWP. In contrast, 
in our approach a systematic determination of all DWP 
within a specific parameter range has been conducted. 
This may hamper the comparison of DWP, obtained by 
the two different methods, even further. 

Having in mind the dominance of the small B-particles 
in the transition between two minima it becomes obvi- 
ous that the first two definitions of the participation ratio 
better reflect the dominance of the single-particle char- 
acter of DWP transitions. Thus we feel it is more in- 
tuitive to speak of 3 rather than of 7-10 particles which 
dominate the translation between DWP relevant for the 
low-temperature anomalies. In any event, from a the- 
oretical point of view it it is the whole distribution in 
Fig. IHl which fully characterizes the nature of the transla- 
tional dynamics between minima in the low-temperature 
regime. 

One might, of course, ask whether this result is spe- 
cific to our binary Lennard- Jones system. Actually, we 
already saw that the DWP with central A-particles in- 
volve more cooperative processes. Our preliminary re- 
sults for Silica indicate that the participation ratios are 
similar to those of the A-particles. In general one might 
speculate that for TLS, which mainly contain very similar 
molecules, the transition is more collective (like in pure 
Silica) and behave like the A-particles in the present case, 
while for TLS consisting of small molecules in a matrix 
the transition is similar to the more localized process, as 
seen for the B particles. 

Now after validating the methods the next step is to 
get a closer insight into the microscopic properties of the 
DWP and to perform a direct comparison with exper- 
iments on the low-temperature anomalies. Work along 
this line will be published elsewhere. 
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